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ABSTRACT 

We determine the slope of the near infrared extinction power law (Aa cx A~") for 
8 regions of the Galaxy between 1-27° and -100°. UKIDSS Galactic Plane Survey 
data are compared, in colour-colour space, with Galactic population synthesis model 
data reddened using a series of power laws and convolved through the UKIDSS filter 
profiles. Monte Carlo simulations allow us to determine the best fit value of a and 
evaluate the uncertainty. All values are consistent with each other giving an average 
extinction power law of a=2.1A'^^Q^^. This is much steeper than most laws previously 
derived in the literature from colour excess ratios, which are typically between 1.6 
and 1.8. We show that this discrepancy is due to an inappropriate choice of filter 
wavelength in conversion from colour excess ratios to a and that effective rather than 
isophotal wavelengths are more appropriate. In addition, curved reddening tracks, 
which depend on spectral type and filter system, should be used instead of straight 
vectors. 

Key words: ISMiHII regions. Starsiformation. Surveys. 



1 INTRODUCTION 

Due to its affect on most astronomical observations, inter- 
stellar extinction remains an important issue for the study 
, of virtually all objects. Stars positioned towards the Galac- 
■ tic Plane, particularly those inside or behind star forming 
regions, are affected most due to the increasing amounts 
of interstellar dust along the line of sight. The extinction 
suffered at visible wavelengths is typically an order of mag- 
\ nitude larger than that suffered in the K band. Furthermore, 
, the wavelength dependence of interstellar extinction varies 
far more in the optical and ultraviolet than in the infrared, 
where it is commonly treated as univ ersal and can be rep- 
resen ted as power law (Aa oc (|Mathisl 11990': 'Draine' 
l2003l ). There are however noted cases in the literature of 
significant variance with flatter extinction laws arising, due 
toCTain growth, in dense clouds (|He et al ] |l995l : lRacca et al.1 
I2OO2I ). 

Due to the lessened effects of extinction, recent Galac- 
tic surveys now use infrared light to probe the densest parts 
of the Galaxy. Although the amount of extinction suffered 
is reduced, extinction correction is still required when at- 
tempting to deredden obscured sources. This is particularly 
true in star forming regions where there can be strong dif- 
ferential r eddenin g betw een n eighbour i ng clu ster members. 

Both iMathid (|l990l ) and IPraind (|2003l ) quote several 
values of a, fitted by several authors, between 1.6 and 1.8. 
These values have mostly been derived by first measuring 
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the near-infrared (NIR) colour excess ratio of well studied 
stars with known spectral types. NIR colour excess ratios 
can be converted to a value of a using the following rela- 
tionship: 
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where Aj is the assumed wavelength of the J photometric 
filter. 

It is commonly noted in the literature that the effec- 
tive filter wavelengths are colour dependent. Using the above 
equation, we would derive different values of a for effective 
wavelengths computed from different colour spectra. How- 
ever the difference between effective filter wavelengths de- 
rived from sources of differing spectral type and reddening 
has usually been considered negligible. With the advance- 
ment of deeper NIR surveys it is now possible to observe 
objects with highly reddened spectra and the previous as- 
sumption must be questioned. 

For highly reddened lines of sight in the Galactic Plane, 
individual stars with well known spectral types are rare. 
Despite this it is still possible to analyse the interstellar 
extinction usin y the red ^ant locus on colour-colour dia- 
grams fCCDsl. ll^ebeto uw et aT (2005") isolated red clump 
giants to determine colour excess ratios and the wavelength 
dependence of interstellar extinction (Aa/Ak) towards two 
different lines of sight. The first is noted as an unremarkable 
region, the second contains the giant HII region RCW 49, 
and a field unrelated to the HII region. They conclude that 
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the wavelength dependence of interstehar extinction is re- 
markably similar desp ite the difference in Galactic location. 
Marshall et al ]i006) created a three dimensional extinction 
distribution model of the inner Galaxy using a catalogue of 
synthetic stars, created using the Be sangon stellar pop ula- 
tion synthesis model of the Galaxy ([Robin et al.) as- 
suming the extinction power law of lMathisI (|l99C)[ ). 

In this paper we derive a NIR extinction power law 
along 8 different lines of sight, distributed along the Galac- 
tic Plane (|b|<5deg), each cen tred on an object i n the 
Red MSX Source (RMS) survey ("Urauhart et"al]l200 /t ) . The 
RMS survey is an ongoing mult i- wavelength observational 
programme designed to return a large, well-selected sam- 
ple of massive young stellar objects (MYSOs) with the pri- 
mary basis from the Midcourse Space Experiment (MSX). 
For this reason, towards the centre of each region of sky 
studied, there will be a massive star forming region, con- 
taining either an MYSO or an ultra-compact HII (UCHII) 
region. Therefore the data will also contain highly reddened 
objects, situated inside or behind the thick natal cloud of 
the massive star. However, since the regions of sky studied 
are large, the dominant cause of extinction will still be due 
to the diffuse ISM. 

Each source in the RMS survey has a computed kine- 
matic distance, from ^^CO observations (l Urauhart et al.l 
I2OO8I ). used to determine the luminosity of each MYSO can- 
didate. This however is not without certain pitfalls. Kine- 
matic distances are typically only accurate to ~lkpc, and 
any RMS object positioned in the i nner Galaxy w il l pos - 
sess a near/far kinematic ambiguity. iBusfield et al.l (|2006l ) 
attempted to solve this ambiguity, to varying degrees of suc- 
cess, using the HI self-absorption technique to determine the 
near or far distance. In order to calculate an independent dis- 
tance to clusters associated with MYSOs, Stead & Hoare (in 
prep.) fit isochrones to the intrinsic colours of individually 
dereddened cluster members. Although the amount of ex- 
tinction suffered at infrared wavelengths is reduced, in com- 
parison to optical wavelengths, extinction correction is still 
required when attempting to deredden obscured sources. 
This is particularly true in star forming regions where there 
can be strong differential reddening between neighbouring 
cluster members. Small changes in the extinction power law 
used to deredden individual objects can have a large effect 
the results. It is therefore necessary to determine an accu- 
rate value of a specific towards a particular line of sight. 

In section [2] we describe the NIR data used to derive 
the extinction power law. In section [3] we discuss the im- 
plications of allowing the filter wavelengths of photometric 
systems to evolve with the progressively reddened observed 
spectra. In section |4] we describe the method used to derive 
a value of a. In section [5] we derive NIR extinction power 
laws for the remaining 7 regions. We draw comparisons be- 
tween our work and that of previous authors in section [6l 
We summarise our conclusions in section [71 



2 DATA 

We extract data for each region from two different NIR 
catalogues, the United Kingdom Infrared Deep Sky Survey 
Galactic Plane Survey (UKIDSS GPS) and 2MASS. 



2.1 The UKIDSS GPS 



The U KIDSS Galactic Plane Survey (GPS, see |LucaselaD 
(|2008l )) covers the region of the Galactic Plane accessi- 
ble by the United Kingdom Infrared Telescope (UKIRT), 
in t he J(1.248Mm ), H(1.631/im) and K(2.201/im) fil- 
ters (|Lawrence et al.1 l200/t ). The GPS data are obtained 
from the Wide Field Camera (WFCAM) Science Archive 
(|Hamblv et al.l lioosl ) . The median 5- cr depths in the sec - 
ond data release of the GPS, given bv lWarren et al.1 (|2007h . 
are J=19.77, H=19.00 and K=18.05 (Vega system). How- 
ever the survey depth is spatially variable, typically decreas- 
ing longitudinally towards the Galactic centre and latitudi- 
nally towa rds the G alactic Plane as the fields become more 
crowded. Lucas et al. (2008) determine the modal depths in 
uncrowded fields to be J=19.4 to 19.65, H=18.5 to 18.75 and 
K=17.75 to 18.0. 

The WFCAM Science Archive enables several differ- 
ent types of data quality cuts to make it possible to select 
only the most reliable data. These include ellipticity cuts to 
remove blended objects, warnings to remove sources with 
saturated pixels, astrometry cuts to remove sources mis- 
matched between different filters, and finally photometric 
error cut s. Th e data cuts used are the same as those pre- 
sented in ILucas et al.l (|2008l ) . The only alteration we make 
is a photometr ic error cut o f 0.03 mag in each photometric 
filter, whereas iLucas et al.l (|2008! ) apply a 0.03 mag colour 
cut. The use of such data quality cuts mean that no sources 
fainter than approximately J'-^18.0, H~16.8 and K'-^lG.l are 
present in the data. An example CCD created using UKIDSS 
data is presented in Fig[T] The UKIDSS GPS data saturate 
at K<12 and so we use 2MASS data to provide a way to 
analyse the brighter stars in each field. 



2.2 2MASS 

The second catalogue i s the Two M icron All Sky Survey 
(2MASS, see SkrutskiellaD (|2006lll covering the whole 
sky i n the J(1.235/im) H(1.662/im) and Ks(2.159//m) fil- 
ters ( Cohen et alll2003l \). The 2MASS data have median 
3-cr depths of J=17.1, H=16.4 and K5=15.3. Only 2MASS 
sources with photometric errors in each band less than 0.1 
mag have been selected. A larger photometric error cut was 
used on 2MASS data, in comparison to the UKIDSS data, 
as reducing the cut to that of the UKIDSS data (0.03 mag) 
drastically reduces the number of returned sources. The ap- 
plication of this photometric error cut means that no sources 
fainter than approximately J^16.5, H^15.3 and K^14.9 are 
present in the data. An example CCD created using 2MASS 
data is presented in Fig[T] 



3 EFFECTIVE FILTER WAVELENGTHS AND 
THE EFFECT ON COLOUR EXCESS 

As mentioned previously in Section [1] it is commonly noted 
in the literature that effective filter wavelengths are colour 
dependent but the differences have usually been considered 
negligible. To investigate how the filter wavelengths change 
for sources with varying spectral t ype and Ay, spectra from 
the stellar atmosphere models by Castelli Kuruc3 (|2004h 
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Figure 1. An example CCD of UKIDSS data (bottom) and 
2MASS data (top). The 2MASS data have been offset along 
the J-H axis by 2.5. Due to the increased photometric depth 
of the UKIDSS data, the giant branch is more prominent com- 
pared with the 2MASS CCD. There are 6638 UKIDSS sources 
and 2694 2MASS sources extracted from the same region of sky, 
centred on the RMS source G48. 9897-00. 2992, covering an area 
of 60^x6^ The data cuts detailed in the text have been applied. 
The UKIDSS data (and 2MASS data to a lesser extent) contain 
many sources scattered below the bulk of the data (in each case 
highlighted in red). These sources most likely possess an IR excess 
due to emission from circumstellar dust. Such stars are abundant 
in star forming regions. Curved reddening tracks, created with a 
KOIII spectrum and a=2.1 specific to each photometric system, 
are plotted as solid black lines. Straight reddening vectors, given 
the same initial gradient as each corresponding reddening track, 
are overlaid as dashed lines. 



(hereafter referred to as the CK04 stellar models) are red- 
dened using a range of a, with extinction in K (Ak) that 
varies between and 3 mag in 0.005 intervals. The spec- 
tra are then convo lved through the J, H and K UKIDSS 
(iHewett et al.ll2006l l and J, H and K^ 2MASS ('Cohe n et all 
I2OO3I I filter profiles using t he H ST Synphot programme as 
described in lHewett et alJ (|2006l l. 

By convolving stellar spectra through the filter profiles 
of each photometric system, we are able to model how the 
effective wavelengths depend upon source colour. It is this 
effect that causes the reddening "vector" to appear to curve 
as the spectra are progressively reddened. Fig|2t. illustrates 
the difference between each filter by comparing the ratio 
between the effective wavelength of each filter, created with 
the spectra of an AOV CK04 stellar model at zero reddening, 
and the effe ctive wavelength of each filter as Ak increases. 
iGolavl ([T974II defines the effective wavelengths of photomet- 
ric filters as follows: 

_ fXF{X)S{X)dX 

jF{X)S{X)dX- ^ ' 

The significantly differing trends of the J filters, pre- 
sented in Fig|2jL, cause the Ej-h reddening track to reverse 



its direction of curvature between the two photometric sys- 
tems, as shown in Fig|2]3. We use the term "reddening track" 
as opposed to the term "reddening vector" for two reasons. 
The term reddening vector is used to describe the magnitude 
and direction of the straight line along which a generic star 
moves, from its original position on a CCD, when observed 
through a specific amount of extinction. As the changing ef- 
fective wavelengths cause this path to curve it can no longer 
be referred to as a vector. Also the amount of curvature 
this path possesses is dependent upon the spectral type of 
a given source. Throughout this paper the term "reddening 
tracks" is used to describe a group of tracks created with 
the same value of a but differing spectral type. The value of 
a used to create each reddening track is specified when nec- 
essary. Fig|2]3 compares the difference in Ej_i^, as Eh-k 
increases, between the static (filter wavelengths at zero red- 
dening for an AOV CK04 stellar model) and the evolving 
filter wavelengths of AOV, KOIII and MOIII CK04 stellar 
models as each are progressively reddened in Ak=0.005 in- 
tervals. Had the filter wavelengths not evolved with the pro- 
gressively redder spectra, the reddening tracks created by 
each model, in either photometric system, would have re- 
mained at AEj-H=0. 

Using UKIDSS data, with typical errors of ^-^0.04 in 
colour- colour space, it is not possible to resolve the differ- 
ences between progressively redder spectra, or different red- 
dening tracks, up to about E//_k=1-2 (or Av~15, using 
the reddening law discussed in Section |4T]) . It is therefore 
safe to assume static effective wavelengths for objects suf- 
fering from a low to intermediate amount of extinction. For 
heavily reddened objects, such as those found inside or be- 
hind star forming regions, this assumption of static effective 
wavelengths no longer holds. The progressively redder spec- 
tra in Fig|2]3 display greatly differing amounts of curvature, 
and the curvature of each reddening track is not consistent 
between filters. For instance, when convolved through the 
UKIDSS filters, the MOIII spectra posses the most curva- 
ture, however when convolved through the 2MASS filters it 
is the KOIII spectra that possess the most curvature (in the 
opposite direction). 



4 MODELLING THE COLOUR COLOUR 
DIAGRAMS 

Each NIR extinction power law is derived by comparing data 
from UKIDSS and 2MASS, with a catalogue of synthetic 
stars created using the Besancon stellar population synthesis 
model of the Galaxy. 

4.1 The Stellar Population Synthesis Model of 
the Galaxy 

A complete description of the Besancon model inputs can 
be found in iRobin etHI (|2003l l. however it is summarised 
here for completeness. The model contains four populations 
of stars, thin disc, thick disc, spheroid and bulge. Each pop- 
ulation is described by a star formation rate history, an ini- 
tial mass function, an age or age range, a set of evolutionary 
tracks, kinematics, metallicity characteristics, and includes 
a white dwarf population. The extinction is modelled by a 
diffuse thin disc, described in terms of visual magnitudes per 
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Figure 2. a (left): Illustrates how the effective wavelengths, ^J^ff (^^^ lines), ^H^ff (blue lines) and XK^ff/^Ks^ff (black lines) change 
as the spectra of an AOV CK04 model are progressively reddened. UKIDSS effective wavelengths are plotted as solid lines, 2MASS as 
dashed lines, b (right): Shows the reddening tracks of a progressively reddened (with a power law of a=2.1) CK04 AOV model in red, a 
CK04 KOIII model in blue and CK04 MOIII model in black. Tracks created using the UKIDSS filters are shown as solid lines and those 
created using 2MASS filters are shown as dashed lines. Note the different directions of curvature. 



kiloparsec (mag/kpc.) The model also allows for the inser- 
tion of_dlscrete_£louds_^ a specified Ay and distance. 

[Marshall et~all (|2006l ) use 2MASS and Besancon model 
data to model the Galactic interstellar extinction distribu- 
tion in three dimensions to over 64000 lines of sight, each 
separated by 15^ in the inner Galaxy (hereafter referred to 
as the M06 distributions). In order to estimate the extinc- 
tion along our chosen 8 lines of sight, we attempt to replicate 
the M06 distributions by using a set extinction model to de- 
scribe the thin disc (mag/kpc) and insert discrete clouds at 
specified distances. 

The Besancon model data contain the distance, the vi- 
sual extinction, the apparent J, H and K magnitudes (given 
in the Johnson-Cousins system), and the spectral type (and 
therefore the intrinsic colours) of each synthetic star. The 
Besancon intrinsic magnitudes are converted to the UKIDSS 
and 2MASS photometric systems and reddened to the de- 
sired extinction along the CK04 reddening tracks, using sev- 
eral different power laws (values of a). 

To determine the extinction along each redden- 
ing track, the visual extinction (Ay) in the Besancon 
model data is converted to an extinction in the J 
band (Aj) of each photometric system. For UKIDSS 
Aj/Av=0.2833 and for 2MASS A j/Av=0.2 900. These val- 
ues are calculated using the .Cardelli et al.l (|l98 9) extinc- 
tion curve with Rv=3.1 (obtained from the Trilegal web- 
site ^http: / / stev.oapd.inaf.it / cgi-bin/trilegal ) . The value of 
A J for each synthetic star in the Besancon model data is now 
matched with a value of Aj, and therefore the correspond- 
ing values of Ah and Ak, in the desired reddening track. 
Using these values of Aj, Ah and Ar we can compute the 
apparent magnitudes of each synthetic star, in the desired 
photometric system, using a specific value of a. The limiting 



magnitudes of the Besancon model data are made consistent 
with those of the observed field and photometric system. 
This is beneficial since the photometric depths of each region 
can change depending on the chosen Une of sight. As sug- 
gested in the previous section, stars of varying spectral type 
produce reddening tracks with differing amounts of curva- 
ture, particularly those created in the UKIDSS photometric 
system. This difference is noticeable at large amounts of red- 
dening. For this reason, depending on the spectral type of 
each synthetic star, different reddening tracks created with 
CK04 stellar models of equivalent spectral type, are used to 
redden each star. 



4.2 The modelled regions 

To determine accurate values of a we need to extract data 
from large areas to return an adequate number of stars. How- 
ever the stellar population and therefore intrinsic colours 
may change over large areas. We have analysed the distri- 
bution of the spectral types of giants along several lines of 
sight using Besangon model data. From this we have found 
that the distribution of background giants can vary dra- 
matically over a few tenths of a degree in Galactic latitude 
but remains constant over several degrees in Galactic lon- 
gitude (with the exception of regions towards the Galac- 
tic centre). For this reason our chosen areas are typically 
1=60^ X b=6^ in size. As our method requires a large num- 
ber of sources to determine a power law, in regions where 
the source count drops significantly (i.e. towards the outer 
edges of the Galaxy), we extend the area in the longitudi- 
nal direction to 90^x6^ in size. In order to avoid observing 
the Galactic bar, since this is not modelled in the Besancon 
model data, we do not attempt regions within 27° of the 
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Figure 3. A CCD created from Besangon model data that have 
not been converted from the Johnson-Cousins system and have 
therefore been reddened along a reddening vector. The data 
have been extracted from a region centred on the RMS source 
G48.9897-00.2992 (as in FigH]). There are 6075 sources and the 
apparent magnitude range is consistent with that of the UKIDSS 
data presented in FigH] The blue points cover stars with a spec- 
tral type range KO-KIIII, the red K4-K6III and the black cov- 
ers those with spectral type K8III and later. The giant branch 
appears to show some curvature due to the dominant observed 
spectral type varying with depth through the galactic plane. 



Galactic centre (iGarzon et al.lll997f ). Also, as there are no 
M06 distributions available with a Galactic longitude >100°, 
we only consider regions between 27° and 100° in Galactic 
longitude. 

Fig |3] contains Besancon model data extracted from the 
same region as the data in Fig[T] As the dominant observed 
spectral type varies with depth through the galactic plane, 
so will the dominant intrinsic colours. This causes the giant 
locus to show some curvature even when reddened along a 
reddening vector. The separation of giant intrinsic colours is 
more pronounced in the UKIDSS photometric system than 
in the 2MASS photometric system. This causes the giant 
locus, as apparent in Fig[T] to show more curvature when 
created using UKIDSS data than with 2MASS data. 

4.3 Comparison of model and observed CCDs 

The Besancon model data are reddened using a range of 
power laws and compared to real data in colour-colour space. 
To determine which power law creates the best fit between 
the synthetic and real data we aim to minimise the following 
quantity: 

1 ^ 

^ = - ^^^{{S{H-K)i-R(H-K),)'^ + {Si^j-H),-R(J-H)i)'^)^ , (3) 

where L is defined as the mean separation, in colour-colour 
space, between each real star (R) and its nearest synthetic 
star neighbour (S). We perform Monte Carlo simulations to 



calculate the la deviation of L and from this a value of a, 
with la errors, is determined. We also sample from a real- 
istic distribution of photometric errors in the simulation to 
ensure we are treating the model data in as similar a way as 
possible to the real data. Typical observational photomet- 
ric errors are created using an exponential function, that 
describes the relationship between apparent magnitude and 
photometric error, (e.g. Jerr = A + exp(B x Jmag - C), 
where A, B and C are coefficients specific to each photomet- 
ric system) and a random Gaussian distribution. Unlike the 
real data, the synthetic data do not contain sources with an 
IR excess. To avoid any contamination in the fitting proce- 
dure, potential IR excess sources are removed from the real 
sample before a fit is made. This process is done by eye by 
removing any sources that extend away from the bulk of the 
data, as illustrated in Fig[T] In each of the 8 regions such 
stars account for ^1% of total sources and therefore have 
little effect on the fit. 

As an example we use data, extracted from around the 
RMS source G48.9897-00.2992, shown in FigH] G48.9897 
has been identified as a HII region with a derived kinematic 
distance of 5.2kpc. Initially we first derive extinction laws 
using Besancon model data created with a simple extinc- 
tion distribution models with a fixed amount of extinction 
per distance (mag/kpc). To determine the deviation of L, 
we compare 200 different sets of Besancon model data with 
the real data from the surrounding area. So we require 200 
times as many synthetic stars as there are real stars (in 
this particular case there are ^6,640 UKIDSS sources and 
so '-^1,330,000 synthetic stars are required) to draw 200 dif- 
ferent random selections of stars, each containing a similar 
number of stars as the real Galactic field, in order to com- 
pare with the real data. A large amount of Besancon model 
data are extracted from a pencil beam centred on the RMS 
source, thereby effectively giving every synthetic star the 
same 1 and b as the RMS source. As mentioned in Section 
14.21 we avoid the problem of the changing stellar population 
by restricting the range of Galactic latitudes used to extract 
the UKIDSS and 2MASS data. 

FigU^ illustrates how L varies with a using one ran- 
domly drawn set of Besancon model data created with an 
extinction distribution model of 1.5mag/kpc. The smallest 
value of L occurs when the Besancon model data best fit 
the real data, and in this particular case this occurs when 
a=2.150. This process is repeated 200 times to determine 
the distribution of L for each value of a. Fig|4]3 contains a 
histogram illustrating this distribution, with a skewed Gaus- 
sian fitted to determine the median and la deviation. The 
value of a that creates an L distribution with the smallest 
median defines the final value of a. The random error is de- 
fined by the values of a where the la deviation subtracted 
from the median value of L, is smaller than the smallest me- 
dian. This process is illustrated in Fig|4]: where the values of 
a=2.100 to 2.175 are within a la deviation of the smallest 
median. 

The fitting process is illustrated using Besancon model 
data, reddened using three different extinction laws along 
the CK04 reddening tracks, and UKIDSS data in Fig|5l 
The final value for UKIDSS data, using an extinction dis- 
tribution model of 1.5mag/kpc, is a=2.150ibo;o84 5 these er- 
rors include accounting for the discretisation of a. The pro- 
cess is repeated using 2MASS data deriving a final value of 
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Figure 4. a (left): a vs L using one randomly drawn set of Besangon model data. The smallest value of L occurs when the Besangon 
model data best fit the real data, b (centre): The distribution of L for 200 different sets of Besangon model data when a=2.15. Also 
plotted are the fitted Gaussian and the median, c (right): The values of L where the Icr deviation error bar pass below the dashed blue 
line (smallest median) define the random error in a. 





Figure 5. a (left): Contains UKIDSS data (red points) plotted against Besangon model data (blue points), reddened along the CK04 
reddening tracks created with q;=2.00. The UKIDSS and Besangon model data are presented twice with an offset along the J-H axis and 
the overlay switched, b (centre): As (a) but with a=2.15. c (right): As (a) but with a=2.30. Unlike the UKIDSS data, the Besangon 
model data do not contain any sources with an IR excess. 



a=2.025±o;o8i • These two values of a agree within la un- 
certainties and give an average law of a=2.088=bo;o59- This 
value of a is significantly larger than the majority of val- 
ues present in the literature. All derived extinction laws are 
presented in Table [1] 



4.4 The extinction distribution 

To test how dependent the derived extinction power law 
is on the extinction distribution model used to create the 
Besangon model data, we perform two tests. First we vary 
the simple extinction distribution model of 1.5mag/kpc over 
a range of intervals between 1.2mag/kpc and 1.8mag/kpc 
corresponding to a change of ±20%. Fig|6^ contains a con- 
tour plot illustrating the simultaneous minimisation of L as 
the extinction model (mag/kpc) and value of a vary. The 
elongated shape in the direction of the y-axis suggests that 
the fitting technique is not strongly affected by the change 
in extinction per kpc. Table [1] contains the laws derived us- 
ing the 1.2mag/kpc and l.Smag/kpc extinction distribution 
models, in both photometric systems. Each is consistent to 



within Icr errors of the laws derived using the 1.5mag/kpc 
extinction distribution model. 

For the second test we randomise the amount of ex- 
tinction each source suffers in order to represent the patch- 
iness of extinction over the observed region. The maximum 
amount of extinction suffered by any source in the Besangon 
model data for the region G48.9897, reddened with an ex- 
tinction model of 1.5mag/kpc, is Av-^20. An amount of 
extinction, randomly selected from a Gaussian distribution 
with mean Av=0 and cr= lOmag, was added to (or sub- 
tracted from) each star. In such cases where the subtraction 
of a random amount of extinction would result in Av<0, 
the original value of Ay in the Besangon model data was 
used. This test has produced some realistic features of a 
CCD and some unrealistic features. The CCD now contains 
a small number of sources that have extended away from 
the bulk of the data with H-K>1.8 (FiglBJ)) similar to the 
real UKIDSS data. However, now there is also an over abun- 
dance of sources with a very low Ay. This creates the blue 
spike in Fig|6]3 that appears between H-K '-^-O.l and ~0.1. 
Although this feature is unrealistic, the small amount of ex- 
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Alpha 




Figure 6. a (left): Simultaneous minimisation of L as the extinction model (mag/kpc) and value of a are varied, b (right): Contains 
UKIDSS data (red points) plotted against Besangon model data (blue points) created by randomising the visual extinction each source 
suffers (see text) and reddened along the CK04 reddening tracks created with a=2.225. 






Figure 7. a (right): A distance- Ay plot containing the extinction distribution model, lying within the constraints of the M06 data, that 
is used to derive each extinction power law for the region G48.9897. b (centre): A CCD created using UKIDSS (red points) and Besangon 
data (blue points) reddened with the derived law, a=2.15. c (left): As (c) but with 2MASS data and a=2.05. 



tinction suffered by these sources means their presence will 
have little effect on the fit. Due to the scatter created by 
these randomised data, the error in the derived extinction 
power law is noticeably larger than the law derived using an 
extinction model of 1.5mag/kpc and UKIDSS data. Due to 
this extra scatter, it was not possible to constrain a value 
using 2MASS data. Table [T] contains the derived extinction 
power law and it is consistent with the extinction power law 
derived using the 1.5mag/kpc extinction model. 

For the region G48.9897, the UKIDSS and 2MASS data 
generate CCDs that appear smooth and somewhat uniform, 
being that there do not appear to be any large clumps or 
gaps in the data. It is as though all sources along the line 
of sight roughly suffer from a constant amount of extinc- 
tion per kpc. This is why the Besancon model data created 
with a simple extinction distribution model, of 1.5mag/kpc, 
produces such a realistic looking CCD. Data extracted from 
many of our other regions often do not produce such smooth 
and uniform CCDs. This is because the molecular cloud sur- 
rounding each star forming region analysed often creates a 



sudden jump in extinction, affecting stars inside and behind 
the cloud producing CCDs with gaps in the data. For this 
reason it is not always possible to assume sources in the 
field suffer from a constant amount of extinction per kpc. 
In such cases a more sophisticated extinction distribution 
model will be needed to detail how Ay varies with distance. 
As discussed in section |4l we use the M06 distributions to 
describe the extinction along the line of sight. We use the 
M06 distribution that is spatially closest to the RMS source 
to attempt to match the detailed morphology of the CCDs 
created using UKIDSS and 2MASS data. Each of the M06 
distributions are replicated by using a set extinction model 
to describe the thin disc (mag/kpc) and discrete clouds in- 
serted at specific distances. 

A distance vs Ay plot for the region G48.9897 is pre- 
sented in Figl^L containing the M06 data and extinction dis- 
tribution model used to derive extinction power laws, consis- 
tent with the laws derived using the 1.5mag/kpc extinction 
model. The new extinction distribution model possesses a 
sharp rise in Ay at '^5kpc; this is consistent with a massive 
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star forming region situated at the RMS kinematic distance 
of 5.2kpc along the hne of sight. There is only a small im- 
provement in the size of L which is likely to be why the 
exact same value of a results in the best fit for the UKIDSS 
data. The size of the error however is reduced slightly. Table 
[T] contains the derived extinction power laws for each pho- 
tometric system. 

The 1.5mag/kpc extinction model is also plotted in 
FiglZt- Despite the constant amount of extinction suffered 
per kpc, the 1.5mag/kpc model (black line) begins to curve 
after ~12kpc. This effect is present (to a greater or lesser 
extent) in the majority of the M06 models used in this pa- 
per. It is an output of the Besancon model and occurs as 
the edge of the synthetic Galaxy, and therefore the end of 
the thin disk, is approached. Fig|7l3-c display the resultant 
CCDs using UKIDSS and 2MASS data, created using the 
final extinction distribution model. Note a slight deficit of 
model data (blue points) in Fig|7l3 around H-K'^1.0, in com- 
parison to FiglSJ), due to the increase in the extinction dis- 
tribution. 

To test the dependence of the results on the assumed 
intrinsic colours, a second set of intrinsic colours and redden- 
in g tracks were c reated from the stellar spectral flux library 
bv lPicklesI ([l998h (hereafter referred to as P98). The process 
has been repeated, using all the extinction distribution mod- 
els detailed in this paper, with the P98 intrinsic colours and 
reddening tracks. The results derived for each of the extinc- 
tion distribution models, using both UKIDSS and 2MASS 
data, were consistent with results derived using the CK04 
intrinsic colours and reddening tracks. 



5 RESULTS 

Results for the remaining 7 regions are presented in Table 
[1] Distance-Ay plots including both the M06 data and any 
extinction distribution models used, along with the resul- 
tant CCDs created using UKIDSS data, are presented in 
Figures [8l to 1101 There is one exception with the final region 
G97.9978. Due to the low amounts of extinction suffered and 
the low number of sources used to determine the fit, caused 
by the galactic position of the RMS source, it was not possi- 
ble to constrain a value of a using 2MASS data and the ex- 
tinction distribution model described by the M06 data. The 
M06 data appear to be an overestimate of the amount of ex- 
tinction suffered, and due to the large area covered (90^x6^) 
it is probable that the M06 data was dominated by the high 
levels of extinction surrounding the RMS source. As the 
majority of data (excluding the highly reddened scattered 
sources) appear smooth and uniform we have performed, in 
the same manor as Section 23] a simultaneous minimisation 
of L as the smooth extinction model (mag/kpc) and value 
of a vary. Instead of a distance- Ay plot, FiglTOk contains 
a contour plot similar to Fig|6]3. The contour lines are not 
as smooth as Fig|6]3 and this is reflected by the size of the 
derived errors. 



6 DISCUSSION 

Creating an average of the derived extinction laws for all 8 
regions results in a value of a—2.1A'^^Q^^. All 8 regions agree 



to this average to within la uncertainties, suggesting that 
the value of a for the ISM is constant over the Galactic lon- 
gitude range 27° < 1 < 100°. Furthermore, there appears to 
be no change with differing Galactic latitude or total visual 
extinction. The slope towards each massive star forming re- 
gion is steep so there is no sign of grain growth causing a 
flattening of the law. However the extinction is still likely 
to be dominated by the diffuse ISM in the large areas we 
analysed. 

The ability to test larger Galactic latitudes using this 
method will be limited, as at larger Galactic latitudes there 
is less total extinction and fewer background stars away from 
the Galactic Plane. As shown with the large error in the de- 
rived values of a for the regions G69.5395 and G97.9976, 
the method suffers from a limited amount of sources and 
extinction. 

Mathii (|199Q| ) and lDraind (l2003l ) present the results of 
several authors, many of whom have calculated values of a 
using colour excess ratio methods typically between 1.6-1.8. 
The main difference with this work is that the derived red- 
dening law presented in this paper has been derived using 
filter wavelengths that evolve with the changing spectra of 
pro gressively reddened objects. 

Moore et all (|2005l ) explore the near-infrared extinction 



law in regions of high Ay towards nine UCHII regions. They 
use hydrogen recombination lines to derive a value of a to- 
ward each region. Such a method does not use broad-band 
filters and is therefore not subject to biases caused by evolv- 
ing effective wavelengths. The results of these nine UCHII 
regions give an average value of a=1.59=b0.30, although 
the range of values covered is large. Excluding results with 
T(Br7) > 3 (which corresponds to Av'^25, an approximate 
upper limit in extinction along the lines of sight presented 
in this paper), and only including distant sources (>3kpc), 
thereby making it likely that the extinction is dominated 
by the diffuse ISM, four sources remain. These data give an 
average value of a=1.86=b0.16. This result is consistent with 
the value derived in this paper at the 2a level. 

When colour excess ratios have been used to derive val- 
ues of a, the value depends on the specified wavelengths 
input into equation [T] It is unclear in the literature how ex- 
actly the filter wavelengths have been derived. Some authors 
have just referenced the photometric system and presented 
approximate wavelengths to specific filters (IWhittetl [19881 : 
Martin k Whittet 1990; Mathis 1990). When defining pho- 
tometric systems, often it is either the isophotal or effective 
wavelengths, derived from specified spectra, that are given. 
Isop hotal wavelen gths are derived from the following equa- 
tion taolavlll974l ): 



jF{X)S{X)dX 
J S{X)dX ' 



(4) 



where F(A) is the flux density of the specific spectra and 
S(A) is the throughput, taking into account the detector 
quantum efficiency, filter transmission, optical efficiency and 
the atmospheric transmission. The isophotal wavelength is 
then defined as the wavelength at which the monochromatic 
brightness (or magnitude) of the source spectrum equates 
to the weighted average flux across the filter band. As such 
it is therefore dependent upon the spectrum of the observed 
source. 
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Figure 8. a (top left): A distance- Ay plot containing the extinction distribution model, lying within the constraints of the M06 data, 
that is used to derive each extinction power law for the region G29.9564. b (top right): A CCD created using UKIDSS data extracted from 
the region G29.9564 (red points) and Besangon data (blue points) reddened with the derived law, a=2.20. c (centre left): A distance-Ay 
plot containing the extinction distribution model, lying within the constraints of the M06 data, that is used to derive each extinction 
power law for the region G42.1268. d (centre right): A CCD created using UKIDSS data extracted from the region G42.1268 (red points) 
and Besangon data (blue points) reddened with the derived law, a=2.30. e (bottom left): A distance-Ay plot containing the extinction 
distribution model, lying within the constraints of the M06 data, that is used to derive each extinction power law for the region G53.6185. 
f (bottom right): A CCD created using UKIDSS data extracted from the region G53.6185 (red points) and Besangon data (blue points) 
reddened with the derived law, a=2.175. 
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Figure 9. a (top left): A distance- Ay plot containing the extinction distribution model, lying within the constraints of the M06 data, 
that is used to derive each extinction power law for the region G69.5395. b (top right): A CCD created using UKIDSS data extracted from 
the region G69.5395 (red points) and Besangon data (blue points) reddened with the derived law, a=2.25. c (centre left): A distance- Ay 
plot containing the extinction distribution model, lying within the constraints of the M06 data, that is used to derive each extinction 
power law for the region G77.9637. d (centre right): A CCD created using UKIDSS data extracted from the region G77.9637 (red points) 
and Besangon data (blue points) reddened with the derived law, a=2.20. e (bottom left): A distance- Ay plot containing the extinction 
distribution model, lying within the constraints of the M06 data, that is used to derive each extinction power law for the region G82.1735. 
f (bottom right): A CCD created using UKIDSS data extracted from the region G82.1735 (red points) and Besangon data (blue points) 
reddened with the derived law, a=2.125. 
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TABLE n 

DERIVED POWER LAWS FOR EACH REGION 



Table 1. " see Section 0] for details. The presence of two kinematic distances highlights a kinematic distance ambiguity. 



Region 


Kinematic 
Distance 


Field 
Size (Ixb) 


Extinction 
Model 


Photometric 
System 


No Sources 


Alpha 


Error 


Weighted 
Mean 


G48.9897-00.2992 
(HII region) 


5.2 kpc 


60^x6' 


1.5m/k 

1.2m/k 

1.8m/k 

Randomised 
Extinction" 
M06 


UKIDSS 
2MASS 
UKIDSS 
2MASS 
UKIDSS 
2MASS 
UKIDSS 

UKIDSS 
2MASS 


6,638 
2,694 


2.150 
2.025 
2.200 
2.150 
2.125 
2.025 
2.225 

2.150 
2.050 


+0.060 
-0.085 
+0.181 
-0.081 
+0.059 
-0.109 
+0.129 
-0.204 
+0.085 
— 060 
+0."204 
-0.129 
+0.106 
-0.131 

+ 0.063 
-0.063 
+0.180 
-0.130 


2.088 

+0.071 
^•^'^ -0.116 

2 075 +0-111 

Z.U^O -0.071 

2 100 +°-°^^ 


G29.9564-00.0174 
(HII region) 


6.2/8.5 kpc 


60^x6' 


M06 


UKIDSS 
2MASS 


6,244 
3,870 


2.200 
2.200 


+0.060 
-0.085 
+0.154 
-0.204 


2 200 +°-°^^ 
z.zuu _o.iii 


G42.1268-00.6224 
(HII region) 


5.0/7.6 kpc 


60'x6' 


M06 


UKIDSS 
2MASS 


4,421 
3,024 


2.300 
2.075 


+0.035 
-0.110 
+0.179 
-0.154 


r> -, QQ +0.091 

^.188 _o.o95 


G53.6185+00.0376 
(YSO) 


1.7/8.4 kpc 


60'x6' 


M06 


UKIDSS 
2MASS 


6,108 
3,237 


2.175 
2.050 


+0.060 
-0.085 
+0.215 
-0.120 


9 -|-|o +0.112 
z.iio _o.o74 


G69.5395-00.9754 
(HII region) 


1.2/4.7 kpc 


60^x6' 


M06 


UKIDSS 
2MASS 


4,101 
2,774 


2.250 
2.050 


+0.203 
-0.228 
+0.277 
-0.302 


Z.iDU _o.i89 


G77.9637-00.0075 
(HII region) 


4.5 kpc 


90'x6' 


M06 


UKIDSS 
2MASS 


5,405 
2,211 


2.200 
2.150 


+0.083 
-0.108 
+0.204 
-0.154 


-17^ +0.110 
^•-L'-> -0.094 


G82.1735+00.0792 
(YSO) 


8.1 kpc 


90^x6' 


M06 


UKIDSS 
2MASS 


3,505 
1,489 


2.125 
2.275 


+0.106 
-0.131 
+0.178 
-0.229 


2 200 +°-^°^ 

Z.ZUU _o.i32 


G97.9978+01.4688 
(YSO) 


1.1 kpc 


90^x6' 


1.4m/k 


UKIDSS 
2MASS 


4,408 
2,650 


2.050 
1.925 


+0.178 
-0.328 
+0.178 
-0.253 


1 QS« +0.126 
±.iyoo _Q 207 



This problem is highlighted further when attempt- 
ing to_coriver^_a^_gwen^^ excess ratio to a value of 
a. Ilndebetouw et aL I EM) derive a colour excess ratio 
Eif-i<:/Ej-K=0.3 5=b0.02, towards t he re g!:ion referred to 
as 1=284° (note, 'indebe touw et all (|2005l ) only presented 
the derived colour excess ratios in their results and did 
not compute a value of a). They have assumed that all 
of the sources used to calculate the colour excess ra- 
tio were K2III and so quote the isophotal wavelengths 
of the 2MASS filters, convolved with a K2III star, to be 
Aj=1.240/im, Ai/=1.664/im and AK=2.164/im. Substitut- 
ing these wavelengths and their colour excess ratio, that 
is equivalent to Ej-i//Ei/-K=1.86±0.11, into equation [1] 
yields a value of a=1.81±0.20. As errors in the isopho- 
tal wavelengths were not provided they have not been in- 
cluded in the calculation. We have determined the effective 
wavelengths of the 2MASS filters using a CK04 K2III spec- 
trum to be Aj=1.244/im, Ai/=1.651/im and AK=2.159/im. 
Inputting these wavelengths into equation [1] yields a value of 
a=2.05±0.20. This value is now consistent with the average 
value derived in this paper. 

An attempt was made to derive an extinction law to- 
wards the region 1=284° using the method described in this 
paper. However only 2MASS data were available as the 



Galactic location is outside the range of UKIRT. Further- 
more, as the region covered 0.6° in Galactic latitude, thereby 
making it difficult to determine an extinction distribution 
model using M06 data, and the levels of total extinction 
were relatively low, it was not possible to constrain a value of 
a using 2MASS data. However, as Indebetouw et al. (200i) 
have assumed that the selected stars used to calculate the 
colour excess ratio are all of the same spectral type, the 
colour excess ratio can be considered to be an approxima- 
tion to the curved reddening track (essentially a straight line 
fitted through the path of any reddening track with a gra- 
dient equal to the colour excess ratio) . As Ilndebetouw et all 
(|2005l l provided the range, in colour-colour space, over which 
the region was analysed, we have calculated the gradient 
of the straight line fitted to many CK04 K2III reddening 
tracks, created with varying values of a, over the same range 
of data in colour-colour space. Therefore using a K2III red- 
dening track, as opposed to a reddening vector to convert the 
colour excess ratio to a value of a, we obtain a=2.10zb0.20. 
This is consistent with the value we derive. 

As a final test to highlight the problem of choosing the 
correct filter wavelengths, we have replaced the reddening 
tracks, used to derive an extinction law towards the region 
G48. 9897-0. 2992 using the smooth 1.5mag/kpc extinction 
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Figure 10. a (left): Simultaneous minimisation of L as the extinction model (mag/kpc) and value of ol are varied, b (right): A CCD 
created using UKIDSS data extracted from the region G97.9978 (red points) and Besangon data (blue points) reddened with the derived 
law, a=2.05. 



TABLE 2 

Derived laws for G48. 9897-0.2992 using a 
1.5mag/kpc extinction distribution model 



Photometric 


Reddening 


Alpha 


system 


method 




UKIDSS 


Tracks 








1 950+°-^°^ 




^eff 


^•-LJU_Q QgQ 


2MASS 


Tracks 


^•u^^_0.081 




^iso 


1 7^:0+0.156 




^eff 


1 975+°-^2^ 

-^•^'^-0.131 



distribution model, with a reddening vector. Results were 
derived for both photometric systems using, for 2MASS, 
the fla t spectrum isophotal wavelengths from ICohen et alJ 
(|2003') and, for UKIDSS , the isophotal wavelengths from 
[^ku naga Vaccal ([20051), and using effective wavelengths 
derived from a CK04 AOV spectrum. A comparison between 
the results is presented in Table 2. If previous authors have 
in the past used flat spectrum isophotal wavelengths to de- 
rive values of a, then this could explain the reason why the 
value derived in this paper is much larger than those quoted 
in the literature. 

It is of interest to compare ou r value of a=2.14 to grain 
models. Data from [Praine Le3 (1984) show that the ex- 
tinction suffered over NIR wavelengths (1.25/im to 2.20/im) 
is caused predominantly by the carbonaceou s component 
rather than the silicate. iDraine Led (|l984[ ) use a power 
law distribution of grain sizes up to 0.25/im to fit the ex- 
tinction law which resulted in a slope of 1.4 from 1.25/im 
to 2.20/im. |K im et al.l ^mi) 

use maximum entropy to de- 
rive the size distribution that showed an exponential cutoff 
at the upper end and gave a NIR slope of 1.8. If w e exam- 
i ne th e properties of single-sized graphite grains from toraind 
(| 19851 ) the NIR slope increases with increasing size from 1.55 
for grains of radii O.Ol/xm, to 2.04 for 0.07/im and 2.39 for 
O.lO/xm. Hence, a further weighting or extensions to larger 
carbon gra ins may be abl e to m atch the result derived in 
this paper. iBorghesi et al.l (|l985l ) find an IR extinction law 



of ~A ^ for amorphous carbon grains. This suggests that 
NIR extinction is not caused by amorphous carbon grains. 



7 CONCLUSION 

We have determined the slope of the NIR extinction power 
law, using two different photometric systems, along 8 differ- 
ent lines of sight, each with a massive star forming region 
towards the centre. Comparisons between colour-colour di- 
agrams (CCDs) of UKIDSS and 2MASS data were made 
with CCDs of Besancon Galactic model data. The Besancon 
Galactic model, combined with the M06 extinction distribu- 
tion models, reproduce the CCDs from UKIDSS and 2MASS 
data well. This enabled us to derive an extinction power 
law specific to each of the 8 different lines of sight. Each 
of the 8 extinction laws are consistent to within la un- 
certainties and show no variation with Galactic orienta- 
tion or varying total visual extinction. An average extinc- 
tion power law of a=2.14lQ Q5, covers the Galactic lon- 
gitude range 27° < 1 < 100°. This value is considerably 
larger than previous values in the literature. It is neces- 
sary to consider the evolution of appropriate photometric 
effective wavelengths with changing observed spectra, par- 
ticularly in deep fields when applying extinction correction 
techniques. In addition this means that straight reddening 
vectors commonly used in CCDs are in fact curved. The 
reddening tracks used in this paper can be downloaded at 
www.ast.leeds.ac.uk/RMS/ReddeningTracks/ . At the proof 
stage we were made aware o f work similar to our own by 
iFroebrich k del Burgol (|2006h . 
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